Unconventional superconductivity on a topological insulator 
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We study proximity-induced superconductivity on the surface of a topological insulator (TI), focusing on 
unconventional pairing. We find that the excitation spectrum becomes gapless for any spin-triplet pairing, such 
that both subgap bound states and Andreev reflection is strongly suppressed. For spin-singlet pairing, the zero- 
energy surface state in the d^-wave case becomes a Majorana fermion, in contrast to the situation realized in the 
topologically trivial high-r c cuprates. We also study the influence of a Zeeman field on the surface states. Both 
the magnitude and direction of this field is shown to strongly influence the transport properties, in contrast to 
the case without TI. We predict an experimental signature of the Majorana states via conductance spectroscopy. 



PACS numbers: 

Topological insulators represent a new state of matter which 
presently is generating much interest fl~|-|4]. While being in- 
sulating in the bulk due to a charge excitation gap, spin- 
dependent conducting channels are formed at the edges or 
surfaces of such materials. These states form as Kramer pairs 
which are topologically protected, persisting in the presence 
of disorder as long as time-reversal symmetry is preserved. 
The allure of topological insulators stems not only from their 
obvious interest from a fundamental physics point of view, but 
also because they may find use in spintronics [5 1. Recent ex- 
periments have observed the surface Dirac states characteristic 
for 3D topological insulators ll6l . 

Another motivation for studying topological insulators is 
that they provide an arena for excitations that satisfy non- 
Abelian statistics: so-called Majorana fermions Q. Elemen- 
tary excitations with non-Abelian statistics form a centerpiece 
in recent proposals for topological quantum computation (8). 
Majorana fermions have been shown to exist as surface states 
at the junction between a superconductor (S) and ferromag- 
netic insulator (FI) deposited on a topological insulator due to 
the proximity effect ll9 fflD . The formation of Majorana states 
at the interface between a superconductor and a topological in- 
sulator lfl2l FPU and at the boundary of superfluid 3 He-B lfl4l 
has also been proposed. Very recently, an experimental study 
reported the observation of doping-induced superconductivity 
in the topological insulator (TI) Bi2Se3 1151 . 

By depositing superconducting materials with an uncon- 
ventional pairing symmetry on top of a TI, an exciting 
prospect of an interplay between the internal phase of the su- 
perconducting order parameter A and Majorana states opens 
up. In this Letter, we investigate how spin-triplet and spin- 
singlet of-wave pairing interact with the environment of a TI. 
We find that spin-triplet pairing universally gives rise to gap- 
less excitations and that both bound-states and Andreev reflec- 
tion are strongly suppressed. For spin-singlet pairing, we find 
that the zero-energy surface states in the d^-wave case are 
now Majorana fermions in contrast to the case of the topolog- 
ically trivial case of the high-r c cuprates. Several works have 
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FIG. 1: (Color online) We consider a TI where superconductivity 
and/or magnetic correlations are induced on the surface via the prox- 
imity effect to host materials with the desired properties. 

previously investigated zero-energy vortex core states in the 
context of 3 He and cold atoms |fl6l . While these Majorana 
state are localized in the vortex core, the present Majorana 
fermions spread along the interface to the superconducting re- 
gion. Moreover, we show how the traditional zero-bias con- 
ductance peak (ZBCP) serving as a hallmark of the <4 v -wave 
pairing state lfT9l is substantially modified in the presence of 
a Zeeman field. In fact, the characteristic conductance spec- 
tra of s-wave and li-wave pairing may be interchanged in the 
presence of a time-reversal symmetry breaking field. Also, 
we find that the conductance exhibits qualitatively very dis- 
tinct behavior with respect to the orientation of the field, in 
complete contrast to the topologically trivial case. 

We will employ a Bogolioubov-de Gennes approach to 
obtain the bound-states and transport properties of the sys- 
tem under consideration. Using a Nambu basis ^ — 
(\|/|,\|/|,\^|,\|/|), the general Hamiltonian for the surface of 
a TI reads: 

ft _f Ho(k) A(k) \ 

" V-A*(-k) -Ho*(-k))> {L) 

where //o(k) = Vp(c x k x + O y k y ) — /j (^ denotes a 2 x 2 ma- 
trix). The gap matrix A(k) depends on both the orbital- and 
spin-symmetry of the Cooper pair. For a spin-singlet symme- 
try such as i-wave or of-wave, one finds that A(k) = A(k)io v . 
Diagonalization of Eq. ([T]i then yields the standard eigenval- 
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FIG. 2: (Color online) Plot of the dispersion for the bound-state en- 
ergies for several values of the Zeeman field m z . As m, — > 0, one 
obtains |e| — > A in the .s-wave case and |e| — > in the rf^y-wave case. 
We have set jUs/Ao = 100 ^> /jj^/Af). The energy dispersion changes 
sign when m z changes sign. 

ues e = r\ v / (v F |k|-|V) 2 + |A(k)| 2 , where r\ = ±1, |3 = ±1. 
Turning to the spin-triplet case, where A(k) = (dk • o)ia v , a 
surprising result appears. Diagonalizing Eq. ([T| now yields 
the eigenvalues 

e = TiVf|k|-p v /^ + |A(k)|2. (2) 

when we for concreteness consider a triplet state dk = A(k)z. 
Remarkably, the superconducting order parameter simply 
renormalizes the chemical potential and the excitations re- 
main gapless. Several anomalous properties can be derived 
from the dispersion Eq. Q. By evaluating the corresponding 
wavefunction, one may conclude that Andreev reflection is 
strongly suppressed at the interface to a non-superconducting 
region since there is no gap in the charge excitation spectrum 
that can retroreflect a hole quasiparticle. We have checked that 
for any triplet symmetry the anomalous dispersion Eq. Q is 
obtained. Thus, the results for singlet and triplet pairing dif- 
fer qualitatively in a fundamental way, as the excitations are 
gapped in the former case whereas they remain ungapped in 
the latter case. One may speculate whether a superconducting 
state is supported at all in the spin-triplet case due to the appar- 
ent lack of a gap which offers a net condensation energy. The 
unusual behavior of the spin-triplet symmetry appears to be 
a direct result of the band-structure in the TI, where the spin 
couples directly to momentum through the term a ■ k in the 
Hamiltonian. Since the spin will be parallell to the momen- 
tum, one may note that pairing between equal spins (triplet 
pairing) at k and — k is not possible. Note that this is distinct 
from the case of graphene, where the operator O does not rep- 
resent physical spin, but rather a pseudospin index related to 
the sublattices lfT71[T8l . 

In order to investigate how Andreev reflection is influ- 
enced by unconventional pairing in the environment of a 
TI, we turn our attention to spin-singlet pairing and first 
consider the simplest experimental hybrid structure that can 
probe this phenomenon: a normal metal | superconductor 
(N|S) junction. The scattering coefficient for Andreev re- 
flection is obtained via setting up wavefunctions and con- 
necting them at the interface. In the N region, we 



have y N - [l,e i9 ,0,0] + r c [l,~ e - ie ,0,0] + r h [0,0, l,-e- ie ] 
at x = 0, while y s °c f f ,[e i P,e i (P+ 9 '), _ e i (e / -Y+) je -*r+] + 
t h [l,-e- ie ' ,^- e '-y~\e i ^-^]. Note that we have taken 
into account the possibility of anisotropic pairing such 
as d- wave, and consequently defined e'P = u± — 

y/^l±y/e^\M&W/E), ande^ = A(6 ± )/|A(e ± )|, 9+ = 
0', 8_ = Jt — 8'. A difference in doping level between the N 
and S regions is accounted for by /jn sin 8 = /j^sinS', since in 
an experimental situation the S region is often heavily doped 
(/Js 3> /jaO, in which case one may set 8' = 0. In this case, we 
recover the bound-state solution e = for d^-wave pairing 
which is manifested as a ZBCP flT9l . 

We now provide an argument for why the zero-energy 
bound-state appearing in the d xy -wave case is a Majo- 
rana fermion, in contrast to the zero-energy states real- 
ized in the topologically trivial high-T £ cuprates. The cru- 
cial factor here is the spin-degeneracy of the Fermi sur- 
face in the latter case, whereas for a TI this degener- 
acy is lifted. In both cases, the 4x4 BdG Hamiltonian 
H satisfies a particle-hole symmetry &H(k)& = ~H*(—k), 

where = ( j ~\ |22|. From this property, one may 

prove that if \|/ e = [Hi(k),W2(k),Vi(k),V2(k)] is an eigen- 
function for the eigenvalue e, then ©\|/ e (— k)* = \|/_ E (k) = 
[v*(— k),V2(— k),K*(— k),«2( — k)] is an eigenfunction for 
(— e). For a zero-energy bound state e = 0, one must have 
\|/ £ = \|/_ e , leading to internal symmetry relations between the 
coherence factors such as Mi(k) = v*(— k). The Bogoliubov 
quasiparticle creation operator for this state is constructed 
in the usual way as "{'(k) = Wi(k)cl(k) + K2(k)ct(k) + 
vi(k)c-|-(— k) +V2(k)ci (— k). Thus, we see that the Majorana 
criterium y(k) = y^(— k) is satisfied. Now, the distinction be- 
tween the zero-energy state in the cuprates and the present 
context of a TI is precisely the spin-degeneracy which al- 
lows one to split up the 4x4 BdG equations to two sepa- 
rate 2x2 equations per spin. Due to the band-structure on 
the surface of a TI, the e = solution is not spin-degenerate 
and we obtain only one zero-energy mode. As pointed out in 
Ref. Il22l . this guarantees the Majorana nature of the fermion. 
We reemphasize that this is different from topologically triv- 
ial Nldjcy-wave junctions, where the zero-energy solutions are 
spin-degenerate, i.e. "double Majorana" modes. 

Recent work has demonstrated how Majorana bound-states 
are induced in the presence of a Zeeman-field when contacted 
to a s-wave superconductor lOOIIll- We now wish to investi- 
gate this phenomenon when the superconducting order param- 
eter is unconventional, i.e. we consider a N|FI|<f-wave junc- 
tion. In the c/„-wave case, the spin-singlet order parameter 
reads A(8) = Aocos(20 — %/2)ia y which normally supports 
zero-energy states. Setting up the scattering wavefunctions 
and utilizing appropriate boundary conditions, one is able to 
extract the reflection and transmission coefficients. The chem- 
ical potential /j$ = fi^ = /j is assumed to satisfy /j ^> A in or- 
der to accomodate proximity-induced superconductivity, ex- 
cept in the FI region where /jpj = 0. Proper gating of the 
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different regions depicted in Fig. [T] allow for control over 
the local chemical potential. Let us now define the quan- 
tities V = VF(k y — K)/m z and K = yf (v F k y ) 2 +ml/v F . We 
then arrive at the following general expression for the bound- 
state energy in the limit of vanishing normal-state conduc- 
tance On, i.e. L — > °° where L is the width of the FI re- 
gion: e = |A(9)|sgn{<:}/VIT<?, C = tan[ln(£j3_/J3+)/2i], 
Sl± = sin(29 + 5) + sin(8) ± [sin(25 + 9) + sin(9)]. Here, we 
have defined 8 = — iln(v/i) and £ = — 1 in the s- and d x 2_ Y 2- 
wave case while £ = +1 in the d™-wave case. In the s-wave 
case, this expression agrees with the very recent finding of 
Ref. [EQ. In the dyy-wave case, a zero-energy solution ex- 
ists in the absence of magnetization m z , indicating the pres- 
ence of midgap Andreev bound states. To explore how the 
magnetization influences the bound-state dispersion, we plot 
in Fig. [2] the bound-state energy for both s-wave and d xy - 
wave pairing for several choices of m z . In the .s-wave case, 
the bound-states have a dispersion only near 9 = when the 
Zeeman field is small, \m z \ <C /J. A zero-energy solution is 
seen to be allowed for normal incidence as long as m z is fi- 
nite. In the dxy-wave case, the dispersion is also very small 
for \m z \ <C fj, but in this case it lies almost at e = 0. In- 
creasing m z (m z > 0) in the li-wave case has the important 
effect of accomodating finite-energy bound-states when mov- 
ing away from normal incidence. It should also be noted that 
the chirality of the bound-states are determined by sgn{m z } 
in both the s-wave and of-wave case. To see this, note that 
m z — > (—m z ) leads to 8 — > 8 + Jt. Due to the symmetry rela- 
tions A± (8 + Jt) = -A=p(8) and ]n(A_/A+) = -ln(A + /A_), it 
follows from the definition of C that sgn{£?} oe sgn{m z }. 

The experimental signature of Majorana fermions mani- 
fested as bound-states would be a characteristic behavior of 
the tunneling conductance, as we now show. The normal- 
ized conductance G/Gq may be evaluated by defining G — 

/"0 2 d8cos6[l + |o,(9)| 2 - |r e (9)| 2 ], where r e and r h denote 
the normal and Andreev reflection scattering coefficients, re- 
spectively, while we choose Go = G(|eV| ^> Aq) as usually 
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FIG. 3: (Color online) Plot of the tunneling conductance G/Gq for 
an N|FI|S junction in the .s-wave and d AV -wave case. We have set 
\xLjvp = 1. 



done in experiments. In Fig. [3] we plot the conductance for 
an N | FI | S junction and compare the s-wave case against the 
dxy-wave case. The latter is normally expected to produce the 
well-known ZBCP due to midgap resonant states |[T9ll . One of 
the main aims of this work is to investigate if this hallmark of 
the d xy -wa\e state survives in the present case of a TI. 

As seen in Fig. |3ja), the conductance for the s-wave case 
displays two coherence peaks at e = Ao as usual when m z — > 0. 
Upon increasing m z , the appearance of zero-energy states are 
manifested by a large enhancement of the zero-bias conduc- 
tance. Therefore, the two finite-energy peaks are merged 
into one zero-energy resonance when the Zeeman field m z in- 
creases. We also note that the d x 2_ y 2-wave case is qualitatively 
similar to the s-wave case in Fig. (3ja). Consider now the 
dxy-wave case in (b), where a zero-bias peak is present when 
m z — > 0, in agreement with our previous analytical finding. 
However, the evolution of the conductance spectra are now 
opposite to the s-wave case upon increasing m z : the zero-bias 
peak is split into two finite-energy resonances. In effect, this 
means that the characteristic features in the conductance spec- 
tra of i-wave and cZ-wave superconductors can be completely 
reversed by introducing a Zeeman field in the TI. 

Due to the coupling between spin and momentum in the 
band structure of the surface of a TI, it is interesting to check 
whether the direction of the magnetization influences the con- 
ductance spectra. In a topologically trivial N|FI|<4,.-wave 
junction, one can prove analytically that the conductance is 
invariant with respect to the direction of the magnetization m 
of the FI layer. Increasing the exchange field in the FI region, 
the ZBCP splits in the conventional case [23 1, similarly to Fig. 
|3jb). We here show that in complete contrast to the topo- 
logically trivial case, the conductance now features a strong 
dependence on the magnetization orientation. We consider a 
magnetization in the x- and y-direction in Fig. |3jc) and (d), 
respectively. It is seen that depending on the magnetization 
orientation, the conductance features three qualitatively dif- 
ferent types of behavior. For m || x, G/Gq is invariant upon 
increasing m x . For m || y, the ZBCP vanishes upon increasing 
m y . For m || % the ZBCP is split upon increasing m z . This 
strong sensitivity to the direction of m is a new feature com- 
pared the topologically trivial case which pertains directly to 
the anomalous band-structure of the TI. It may be understood 
by noting that m y shifts the Fermi surface while m z opens the 
energy gap in the FI region. This places strong restrictions 
on how the wavefunction in the S region connects to the FT. 
For sufficiently large m y , the Fermi surface is shifted in such 
a fashion that there are no angles of incidence where surface- 
bound states can be formed any more. We note that an inclu- 
sion of the orbital effect due to the vector potential A simply 
would add a component to the magnetization vector as a result 
of the linear energy-momentum dispersion. 

The predicted results in this work can be tested experimen- 
tally by fabricating a hybrid structure such as the one shown 
in Fig. [T] In terms of actual materials, EuO or EuS might be 
suitable as ferromagnetic insulators in this context lEOl . For 
the of-wave superconductor, a high-T c cuprate such as YBCO 
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would be appropriate. Due to the lattice mismatch between 
the host proximity materials and the TI, the induced super- 
conducting order parameter Ao can be expected to be substan- 
tially reduced in magnitude on the surface of the TI, typically 
in the range 0.1 — 1 meV. Concerning the length of the sam- 
ple, we have set pL/vf = 1. Estimating Vf ~ 5.0 x 10 5 m/s 
and fj — 80 meV ETIl , this corresponds to L ~ 40 nm which 
should be experimentally feasible. 

In summary, we have considered the interplay between 
magnetic order and unconventional superconducting pairing 
on the surface of a topological insulator. We find that the 
charge excitation spectrum is rendered gapless for any spin- 
triplet state, such that both bound-states and Andreev reflec- 
tion are strongly suppressed. For spin-singlet pairing, we find 
that the zero-energy surface states in the cL,-wave case are 
now Majorana fermions, in contrast to the case of the topolog- 
ical^ trivial high-r c cuprates. We have studied how Andreev- 
bound states and Majorana fermions are influenced by the in- 
ternal phase of the superconducting order parameter, and find 
that the ZBCP being the hallmark of the d™-wave state is qual- 
itatively strongly modified in the present context. In particu- 
lar, it is highly sensitive to the magnetization orientation, in 
contrast to the topologically trivial case. Our findings can be 
directly tested through tunneling spectroscopy measurements, 
and we have estimated the magnitude of the necessary exper- 
imental quantities. 
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